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A vibronic-exciton model is applied to investigate the mechanism of enhancement of coherent 
oscillations due to mixing of electronic and nuclear degrees of freedom recently proposed as the origin 
of the long-lived oscillations in the 2D spectra of the Fenna-Mattews-Olson (FMO) pigment protein 
complex [Christensson et al. J. Phys. Chem. B 116 (2012) 7449). Here we reduce the problem to 
a model bacteriochlorophyll (BChl) dimer to elucidate the role of resonance coupling, site energies, 
vibrational frequency and energy disorder in the enhancement of vibronic-exciton coherences as 
well as ground-state vibrational wave-packets, and to identify regimes where this enhancement is 
significant. For an heterodimer representing the two relatively strongly coupled BChls 3 and 4 of the 
FMO complex, we find that the initial amplitude of the vibronic-exciton and vibrational coherences 
are enhanced by up to 15 and 5 times, respectively, compared to the vibrational coherences in the 
isolated monomer. This maximum initial amplitude enhancement occurs when there is a resonance 
between the electronic energy gap and the frequency of the vibrational mode. The bandwidth of 
this enhancement is found to be about 100 cm -1 for both mechanisms. The excitonic mixing of 
electronic and vibrational degrees of freedom leads to additional dephasing relative to the ground- 
state vibrational coherences. Here we evaluate the dephasing dynamics by solving the quantum 
master equation in Markovian approximation and observe a strong dependence of the life-time 
enhancement on the mode frequency. It is found that long-lived vibronic-exciton coherences are only 
" generated when the frequency of the mode is in the vicinity of the electronic resonance. Although 

the vibronic-exciton coherences exhibit a larger initial amplitude compared to the ground-state 
vibrational coherences, we conclude that both type of coherences have a similar magnitude at longer 
population time for the present model. The ability to distinguish between vibronic-exciton and 
ground-state vibrational coherences in the general case of a coupled multi-chromophore system is 
discussed. 
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I. INTRODUCTION 



The observation of long-lived oscillations in the two-dimensional spectra of Fenna-Mattews-Olson (FMO) pigmcnt- 
^ — 1 protein complex [l| has renewed the general interest in quantum effects in light-harvesting and other biological 
. systems. Although similar observations were reported previously 0, ||, the interpretation of the oscillations as 
electronic coherence and the suggestion that such dynamical coherence may play a crucial role in achieving the high 
efficiency of electronic energy transfer in photosynthesis |l| have ushered a wave of research on the role of coherence 
■ in excitation energy transfer processes. 
t-H | Numerous theoretical articles have been devoted to understanding and explaining the long-lived oscillations in 
FMO. They have shown that coherent oscillations last for several hundreds of femtoseconds, but under standard 
assumptions about the properties of the protein environment, none of these works have been able to account for the 
, over picosecond dephasing times of the oscillations observed experimentally. The only proposed mechanism which 
predicts electronic coherence with life times >1 ps is the correlation between environmental fluctuations on different 
pigments Experimental confirmation of correlated distributions of pigment energies have been claimed for 

FMO but neither dynamic nor static correlation of the pigment energies have been found in molecular dynamics 
(MD) simulations of the FMO protein environment [12, Furthermore, dynamic correlation, if present, would lead 
to an effective decrease of the system bath coupling strength, and to a corresponding slowdown of the energy transfer 
rates. Such a decrease of the system-bath coupling strength is at odds with works which find an optimal strength 
of the system bath interaction for transport function within the parameter range used by standard theories |13| . KL4| . 
Also, the energy transfer rates obtained by standard theory are in a good agreement with experimental data [15| . and 
the standard "funnel picture" of energy transfer seems therefore to be well supported. 

The assumption that the experimentally observed coherences would be relevant to the biological function of the 
antenna complexes has also been criticized. Such a proposal requires the excitation with coherent superpositions of 
states created by lasers in the laboratory to be, in some way, equivalent to the excitations under in vivo conditions 
[l7j . i.e. by direct sunlight or via transfer from another antenna system. However, serious objections to this view 
have been raised in the literature, arguing that direct excitation by light from the sun (thermal light) does not lead to 
such coherent excitation [l8r-22]. A recent experiment showed that even under coherent excitation of the chlorosome 
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antenna, i.e. with femtosecond laser pulses, no coherent oscillations could be observed (HJ. This implies that any 
coherence induced by the excitation would decay quickly, and that the energy is transferred to the FMO complex in 
an incoherent fashion independently of the excitation conditions. 

The above discussion suggests that the mechanism explaining the oscillations does not need to have a strong impact 
on the energy transfer dynamics or be relevant to the energy transfer efficiency. The problem should therefore be 
studied from the point of view of the experimental signals rather than the (electronic) excited-state dynamics. In 
other words, the origin of the oscillating experimental signal has to be understood before the signal can be translated 
into claims about the energy transfer dynamics. Recently, Christensson et al. [2~ij proposed that the excitonic 
interaction between electronic and vibrational states in FMO serves to create vibronic states (excitonically mixed 
electronic and vibrational states). Such states have a considerable vibrational character, and at the same time have 
an enhanced transition dipole moment due to intensity borrowing from the strong electronic transitions. It was 
also shown that coherent excitation of the vibronic states produces oscillations in the non-linear signal that exhibit 
picosecond dephasing times. Moreover, the concept of vibronic excitons provides a plausible explanation for the 
observation of correlated distribution of site energies The vibronic-exciton states involved in the long-lived 

coherences are to a large extent composed of different vibrational states on the same pigment. This automatically 
leads to a correlation in the fluctuations of the involved transitions even for a random distribution of pigment energies. 
Finally, the vibronic-exciton model has been shown to provide, compared to the purely excitonic model, more realistic 
relaxation rates, as verified against experimental measurements in two cyanobacterial light- harvesting proteins [25]. 

In a recent work, Tiwari et al. |26| used a similar model to show that the mixing of electronic and vibrational 
DOF leads to an enhancement of the excitation of vibrational coherences in the electronic ground state as well, and 
it was also argued that this effect can explain the long-lived oscillations in FMO. In order to distinguish between 
the coherences observed in the excited-state manifold and the ones originating in the ground-state manifold, we will 
assign the term vibrational coherence strictly to the latter ground-state contribution. There, the states are of pure 
vibrational origin although the ability to excite them is enabled by the mixing of vibrational and electronic states in 
the excited-state manifold. When referring to the excited-state coherences with a strong vibrational involvement, we 
will use the terms vibronic, or more specifically vibronic-exciton coherences. 

In this paper, we turn to a model system (an FMO-inspired molecular dimer) in order to systematically investigate 
how the mixing of vibrational and electronic DOF leads to long-lived oscillatory signal in non-linear optical spec- 
troscopy. The paper is organized as follows. In the next section, we explain the mechanism of coherence amplitude 
enhancement on a toy model. In Section [IIII we introduce the model dimer system with the associated notation, and 
we define a measure of the vibrational character of a coherence in Section IIV1 Section [V] introduces the third-order 
non-linear signal resolved in the coherent 2D Fourier transformed spectroscopy and its relation to the dynamics of 
molecular systems, in particular to its coherent dynamics. After determining the basic properties of the enhancement 
for vibronic (excitonically mixed electronic and vibrational) and the ground-state vibrational coherences in Section 
IVI1 we discuss the results in context of FMO and the low frequency vibrational spectrum of BChl-a in Section IVlII 



II. ENHANCEMENT OF COHERENCE AMPLITUDE BY TRANSITION DIPOLE MOMENT 

BORROWING: A TOY MODEL 

To illustrate the mechanism of enhancement of the coherence amplitude that is suggested in Ref. [24j and studied 
in detail in this paper, let us first consider qualitatively a rather trivial example. Let us imagine a system of two 
molecules, one of which has a forbidden transition to the excited state. In the absence of excitonic coupling between 
the excited states of the two molecules, it is not possible to excite a linear combination of the collective excited 
(eigen) states of the dimer (i.e. the coherence observed by non-linear spectroscopy), because one of the transitions is 
forbidden. If we now switch on the interaction between the monomers, it becomes possible to excite the coherence 
between the new eigenstates of the system (see Fig. Q~|^ for illustration) . We assume the transition dipole moment in 
a form 

£ = M e9 |l)(0| + /!..c., (1) 

where |f) = |ei)|ffa) is the allowed collective excited state of the dimer with the monomer 1 in its excited state |ei) and 
monomer 2 in its ground state \g2) and |0) = |c?i ) |ff2) is the collective ground state. Expressed in the new eigenstate 
basis of vectors |a), the transition dipole moment reads 

A= Yl M<x\l))\a)(0\+h.c. (2) 

ct=l',2' 

The new transition dipole moment to the excited eigenstates \a) reads /i Q o = Meg( a |l)- As will become clear in the 
following sections, the amplitude of the contribution of a coherence between two states |a) and \(3) to the 2D spectrum 
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Figure 1: Illustration of the effect of transition dipole moment borrowing and excitonic mixing of exited states in a dimer 
without considering vibrational states (A) and considering a simplified model of vibrational presence (B). In the case (B), a 
heterodimer is considered with an approximate resonance between the one-phonon state of monomer 1 and the zero-phonon 
state of monomer 2. Panel C: Dimer model with a single vibrational mode per monomer. The monomers 1 and 2 can be in 
their respective ground states \g) or excited states |e). In each electronic state, the monomers can occupy any of the vibrational 
levels corresponding to their respective vibrational modes (qi or q2). 



is proportional to the factor A a p = |/2 Q o| 2 |/^o| 2 - In a dimer, we have only states |1') and |2'}, i.e. only one coherence, 
and we obtain 

A V2 , = \» eg \*\(l'\l)f\{2'\l)\ 2 . (3) 
Because |(1'|1}| 2 + |(2'|1}| 2 = 1, the maximum value of ^4i<2' is equal to jl^egl 4 when the mixing is maximum, i.e. 

I(i'|i>l 2 = i 

In a dimer, the coefficient (l'|l) is given by a sine or cosine of the mixing angle = arctan(-^j), where J is the 
resonance coupling and AE is the difference between energies of the interacting states [27j. When J is small, the 
mixing can only be substantial at the resonance, i.e. when AE — > 0. The observation of coherence in 2D spectroscopy 
depends crucially on the ability to excite the two involved excited states simultaneously. In the above case, it is 
enabled by the resonance between excited states. 

In Fig. [TJ3, we present another interesting case. We extend the previous toy model with vibrational states at each 
monomer, and assume both monomers to have allowed transitions to the excited states. In Fig. [T}3, only the first 
vibrationally excited state in the electronically excited state is represented. The one-phonon state of monomer 1 
interacts weakly with the zero-phonon electronically excited state of monomer 2. Because the coupling is small, the 
interaction between the two zero-phonon states of the interacting monomers can be neglected, and the borrowing effect 
occurs only for the collective states |2') and |3') of Fig. QJ3. If, in addition, we assume the transition dipole moments 
to the one-phonon state to be small (i.e. the vibrational mode has a low Huang- Rhys (HR) factor), we obtain the 
situation similar to Fig. [TJ\. At resonance (i.e. at maximum mixing |(2'|e2)| 2 ~ the coherence amplitude Ay%> in 
2D spectrum reads 

M'2> ~ iMeoSolVeoffol 2 ^- ( 4 ) 

This is to be compared with the amplitude of the purely vibrational coherence on a non-interacting monomer which 
go \ and which is small due the small HR factor. The excitonic interaction thus enhances a 
coherence which has partially both the vibrational and electronic characters. 

The situation studied later in this paper is similar to the model illustrated in Fig. [TJ3. We will study coherences 
involving the zero- and one-phonon states of the same vibrational mode on a given monomer. As demonstrated, the 
transition from the electronic and vibrational ground state to the one-phonon state can be enhanced by excitonic 
coupling with the allowed states of other monomers. This was found in Ref. [24| to be the case for FMO, where the 
enhanced coherences where found to be of a prevailingly vibrational character, thus exhibiting significantly prolonged 
life time over the purely electronic coherences. 
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The excited-state mixing between the vibrational states of one monomer and the electronic (zero-phonon) state 
of the other monomers can, however, enhance also the ground-state bleaching signal [26], which was not included in 
our toy model. The life-time of the coherence has not been discussed within the toy model either, because it is not 
relevant to the illustration of the enhancement mechanism. However, it will be of crucial importance throughout the 
rest of the paper. 



III. DIMER HAMILTONIAN: FULL FORMULATION 



In this paper, we concentrate on the relatively strongly coupled dimer of BChls 3 and 4 of the FMO complex 
[28| . which we will later refer to as monomer 1 and monomer 2, respectively. Needless to say, the interaction of the 
other BChls in the complex with this dimer is decisive for its function as a molecular wire, and it influences the 
properties of the oscillations observed in 2D spectra as well. In Ref. [24|, it was demonstrated that the effect of 
long-lived coherence could plausibly originate from the interplay between the vibrational modes local to BChl 3 and 
the resonance interaction between BChls 3 and 4. By reducing FMO to this dimer, we expect to be able to isolate the 
main quantitative contributions to the observed effect, and to simplify the treatment such that the main conditions 
for the mechanism to take place can be identified. The extension of the results to the entire FMO complex is discussed 
in Section EED 

The studied molecular system will be represented by a dimer Hamiltonian, explicitly treating a single intramolecular 
vibration (with frequency luq) on each of the monomers. Unlike in Ref. [24]], we will first consider the formulation 
with an arbitrary number of vibrational-excited levels, which we refer to as a full formulation (full with respect to the 
number of vibrational levels). For practical calculations, we will reduce the number of vibrational levels to a single 
mode on each monomer. This is justified by the low Huang-Rhys factors used in this work. The dimer Hamiltonian 
in the electronic site basis is defined as 

H s = ff« <g> I' 2 ' + E« ® + j 12 {\ ei )\g 2 )( gi \(e 2 \ + h.c.) , (5) 

where lS n > represents the identity operator on the Hilbert space of monomer n, and 
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(pl + (q n - df) 



(6) 



is the Hamiltonian of monomer n. Here p n and q n are the dimensionless momenta and coordinate operators of the 
vibrational mode on site n, respectively, d is the coordinate displacement, and J nm denotes the inter-site coupling. 
The corresponding level model including the vibrational states and the associated notation is presented in Fig. [Tp. 

This system interacts with a bath of protein degrees of freedom (DOF) which is modeled as an infinite number of 
harmonic oscillators characterized by some continuous spectral density, e.g. J(w) = ^~J a g^5(oj — uj a ). The Hamiltonian 
representing the coupling between the bath modes characterized by the Hamiltonian 



and the dimer reads as 
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The first term in Eq. ([5]) corresponds to the linear interaction between the bath and the electronic DOF and can be 

expressed through the spectral density as K^f 1 = tkij a ga [ll], while the second term in Eq. ([5]) describes a bi-linear 
interaction between the intramolecular mode and the bath. Following the work of Refs. [24j and (26|, we will consider 

this interaction to be weak and set = for all n, accordingly. 

Let us now define the notation used throughout this paper to describe the dimer state. As in the standard electronic 
molecular exciton theory, the excited states of the dimer are formed out of the excited states of the monomers. We 
denote the electronic excited states by a diad (nm) where n (to) stands for the state [g or e) of monomer 1 (2) in 
correspondence with Fig. [Tp. The Hilbert space in which the molecular Hamiltonian will be represented can be fully 
described by states \n u m U '), where the first (second) letter describes the first (second) monomer's electronic state, 
and the index v, v' denote their respective vibrational quantum level, so that e.g. 



\9»9w) = \9M 9) )\92)W 2 (9) ), 



(9) 
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\9»e v ,) = \g 1 )\v[ s) }\e 2 }\4 e) }- (10) 
The eigenstates of the Hamiltonian will be written as linear combinations of the states \n v m v i), i.e. 

I a) = ECmJw)- (11) 

The expansion coefficients obtained by diagonalization of the Hamiltonian 

C v m„, = (a\n v m u >) (12) 

will report on the participation of the particular local states in the eigenstates. The eigenstates will be numbered 
with increasing energy e a , a = 1,2, . . . . 



IV. MEASURE OF THE VIBRATIONAL CHARACTER 



To measure the vibrational character of a coherence, we can use the square of the expansion coefficients. We will 
be interested in the coherence between states \eog v i) and \e%g v '). This coherence is a vibrational coherence located 
on monomer 1, and we ignore the vibrational state of the second monomer. Because the ground-state vibrational 
quanta of the second monomer cannot be excited from the thermal ground state by a single interaction with the 
exciting pulse in the optical domain, we can set v 1 = 0. We are interested in identifying which coherences between 
the eigenstates have the character of this local vibrational coherence, and we want to define accordingly a measure 
Xa/3 °f the local vibrational character of the coherence p a g. For a given coherence to be vibrational on monomer 1, 
the state \a) needs to have a character of the zero-phonon state on the monomer 1, i.e. it has to be predominantly 
composed of the state |eo<?o)> an d the state \j3) has to have predominantly the character of the state |ei<7o) ( or w * ce 
versa). This composition is measured by the square of the corresponding expansion coefficients c" m , and m 
which also define the probability of finding the system in the corresponding local states, should we attempt such a 
measurement. The character of the coherence corresponds to the conditional probability of finding the state |eo<7o) by 
measuring on state \a), and of simultaneously finding the state jeiffo) by measuring on state or vice versa. This 
leads us naturally to the definition 

V a = v ( e o9o,eiff ) I a |2| 13 |2 , | a |2 B |2 /io\ 
X.a/3 - X a/ 3 - \ C e g \ C e igo \ +l C ei 3o l l C eoffol ■ 

This quantity has its maximum when the eigenstates are each entirely composed of the zero- or one-phonon states of 
the excited molecule, and thereby correspond to the site basis states (xa/3 = 1 for |a) = \eogo) and |/3) = \e\ga), or 
vice versa). When the states \a) and \/3) are equal mixtures of these local states, so that |Cg ogo | 2 = \c" igo | 2 = |cf o9o 2 — 
l c eis | 2 = 2"' then we have Xa/3 — \- Therefore, coherences p a p with Xap > \ will be considered to have prevailingly 
local vibrational character. We compared the local vibrational character defined by Eq. (|13j) with the composition 
of the eigenstate vectors |c" TO J 2 , and we verified that our conclusions remain the same regardless of the quantity 
used for the definition of the coherence character. Note that the character Xap is a time-independent quantity. It 
refers to eigenstates of the Hamiltonian and their representation in the basis of states local to the chromophores. 
The coherences have life times dependent of the interaction of the system with the bath, but their local vibrational 
character Xa/3 only depends on the system Hamiltonian and does not evolve in time. 



V. 2D NON-LINEAR SPECTROSCOPY AND LIOUVILLE PATHWAYS 



Non-linear spectroscopy has been widely used to study the dynamics of excitonic energy transfer in light-harvesting 
systems, in particular because it is sensitive to the time dependent redistribution of the populations among excited 
states. Details about this class of measurement techniques are outside the scope of this work and can be found 
elsewhere in literature (see e.g. (29l - [3l| ). Two-dimensional (2D) coherent spectroscopy, which is one of the recent 
additions to this class of techniques, has enabled us to directly observe the coherent components of the excited-state 
dynamics [H, 113 • In 2D coherent electronic spectroscopy, the time dependent signal (the detection time is usually 
denoted as t% here) is generated from the interaction of the sample with three consecutive ultrafast laser pulses with 
wave vectors k\, and k$ separated by time intervals t\ and t2 [32l l34l|. The signal is heterodyne detected, i.e. 
the generated field rather than its intensity is measured, and it is spectrally resolved in frequency W3 (£3 and C03 are 
correspondingly related by Fourier transform). The measured field is then numerically Fourier transformed in time 
t\ {t\ -4 wi) resulting in a signal dependent on two frequencies C03 and ui\ and one time delay £2- 2D spectra are 
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Figure 2: Liouville pathways and their respective contribution to the 2D spectrum. (A) Double-sided Feynman diagram describ- 
ing the four Liouville pathways along with their dipole pre-factor, the frequency of the involved coherence and the frequency 
of the exciting and probing pulses. (B) Position of the signal in the 2D spectrum and its characteristics originating from the 
different pathways. Solid symbol denote non-oscillating contributions in t%, open symbols denote oscillatory contributions with 
frequencies u> a p (circle), w gogv (diamond). Signals denoted by circle symbols do not involve any vibrational level whereas those 
denoted by diamonds do (for v > 1). The diads in (B) represent the eigenstates |a), |/3) involved in each signal contribution. 



represented as the frequency- frequency correlation plots at various time delays £2- The 2D signal is related to the 
third-order polarization of the molecular system generated by the exciting laser pulses, and its theoretical description 
is most conveniently based on the non-linear response function formalism in the third-order of perturbation theory 
(see e.g. Ref. [35]). In this section, we will briefly present the various signal components that can be distinguished 
in a 2D spectrum, and we will give some of their characteristics, such as positions in the 2D spectrum, oscillation 
frequency, initial- and time-dependent amplitudes. 

Depending upon the time ordering of the light-matter interactions in the third-order response functions, four 
different types of so-called Liouville pathways (denoted i = 1, . . . , 4) can be distinguished in the generated signal. 
They are represented in Fig. by double-sided Feynman diagrams [35||. Pathways R\ and R4, i.e. with negative 
ti, are of the so-called non-rephasing character, while R2 and R3 are of the rephasing character. Rephasing and 
non-rephasing signals can be measured separately. The Liouville pathways can also be grouped together according to 
the electronic band in which the coherences observed in 2D spectra originate. In R\ and R2, the observed coherences 
stem from the electronically excited state, while for R3 and R4, they come from the ground state. Thus in both 
rephasing and non-rephasing pairs of Liouville pathways, we observe one ground- and one excited-state coherence. It 
is notoriously difficult to distinguish ground- and excited-state-originated signal experimentally, unless one can find 
some secondary characteristics (such as signal life time due to relaxation or frequency) which distinguishes the ground- 
from the excited-state contributions. 

Fig. [2)3 illustrates the position, origin (excited or ground state) and oscillating frequency of the signal generated 
by the different pathways in the 2D spectrum. The non-rephasing pathways _R( ay8,9 ") and R^^' 9 "^ result in peaks 
at the position (w a9D , u agv ) and (uj a g , up go ), respectively, whereas the rephasing pathways i? 2 " /3 ' 9 "' ) and R 1 ^* 13 ' 9 ^ 
generate peaks located in similar positions, i.e. (Lj aga , ujp 9v ). It is interesting to note the role of the index v in the 
pathways. For v = 0, there are no ground-state coherences, while excited-state coherences can appear. The frequency 
of oscillation of the peaks coming from E^ Ptau) and R^ M is uj a @ , while it is given by the ground-state vibrational 
frequency w Sl/So for R^P* 9 ^ anc f R^' 9 "' . For v > 0, the amplitudes of R\ and -R4 pathways are equal, and the same 
holds for i?2 and R3 pathways. Thus, if we would ignore the effect of the coherence life time, there would always be 
the same amplitude of the ground- and excited-state oscillating contributions (y > 0), plus one contribution from the 
excited state (v — 0). 
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Let us now discuss the (initial) amplitudes of oscillating features in 2D spectra. The initial (t 2 = 0) amplitude of 
the third-order response function is proportional to a transition dipole moment pre-factor of the response function 
Ri averaged over an isotropic distribution of orientations of the molecules in the sample. The dipole pre-factors A^' 
for each Liouville pathway Ri are defined in Appendix [X] For the non-rephasing response function R\ , it reads in 
particular 

A al,g» = (OVa ' e i)^09, ' e 3)(/%,/3 • e 2 )(n aga ■ ei)) n (14) 

where e$ represents the unitary orientational vector of the i th laser pulse, lJ. a g„ is the eigenstate transition dipole 
moment and (■ ■ - )q stands for the averaging over an ensemble of randomly oriented complexes. The disorder averaging 
can be accounted for by a constant orientation factor VL a p = i (2 cos 2 9 a p + l) , where 9 a p is the angle between the 
transition dipole moments of eigenstates a and /3 [36] . The amplitude pre-factor then reads 

\n aP . (15) 



Here, \fi\ — y /i 2 + /i 2 + /i 2 is the length of the transition dipole moment vector. 

The time-dependent response function is composed of the dipole pre-factor and the time-dependent part representing 
time evolution of the state of the molecular systems during the time interval between the interactions with light. The 
response functions Ri for all pathways are presented in detail in Appendix^ In particular, the non-rephasing response 
function R\ reads as 

flSL(*3' = ( A ^l^ G ^(h)G' ap (t 2 )G ago (t 1 )) A (16) 

where G agv (i) and G' a p{t) are the evolution operators of the optical coherence p ag „ and of the excited-state coherence 
p a fi, respectively. Here, (. . . )a denotes averaging over the energy disorder. The details of the evolution operators are 
given in the Appendix. The contribution of the corresponding coherence to the signal in the 2D spectra is obtained 
after Fourier transformation of the evolution propagators during the time intervals t\ and t^: 

S^l^M^i) = (<L G ag „(u> 3 )G' af) (t 2 )G ago (u Jl )) A (17) 

with G agv {u>) = Jo +0 ° dt e l ( w ~ w «9,J t - r = t ~7o ! Qg(t) an( j w jj ere Wlj U}3 are ^\ ie excitation and probing frequencies, respec- 
tively. The signal for the non-rephasing pathway i?4, which involves the electronic ground-state coherence can be 
obtained in the same way. Full expressions are given in the Appendix. 

VI. FMO-INSPIRED HETERODIMER IN ONE PARTICLE APPROXIMATION 

In order to quantify the enhancement of coherence amplitude by borrowing of transition dipole moment, as il- 
lustrated in Section UH we study a simplified molecular dimer model. In this section, we will first study the initial 
amplitude of the dipole pre-factor for selected coherences as a function of the molecule inter-site coupling J and the 
energy gap AE. We will include one particular nuclear mode (ljq = 117 cm -1 ) at each of the two sites. Then, we 
will investigate the dynamical behavior of the signal as a function of the dimer energy gap and the vibrational mode 
frequency. Possible consequences of our results for the case of FMO molecular aggregate will be discussed in Section 

rvni 

As mentioned above, the studied dimer is inspired by the FMO lowest energy BChls 3 and 4 which are relatively 
strongly interacting (site coupling Jo = —53.5 cm" 1 ) and form a heterodimer with an energy gap AEq — 110 cm" 1 
(later referred to as the reference energy gap) (28j. The directions of the transition dipole moments were taken from 
the Protein Data Bank file 3ENI j3?J (we used 6*Bchi 3-4 = 107.2°). A Huang-Rhys factor of S = 0.05 has been used. 
In order to assess the influence of the nuclear modes on the dimer spectra, we studied the dependence of the spectra 
on selected model parameters for a single vibrational mode. We study two cases, with frequencies uiq = 117 cm" 1 
and 185 cm -1 . For some parameters and spectral characteristics, we study the system properties as a function of the 
mode frequency in an interval with a width of roughly 200 cm" 1 . 

Only one vibrational level (u max = 1) of a single mode is treated explicitly in the molecular Hamiltonian. This 
is motivated mainly by the small value of the Huang-Rhys factor. We have verified that we obtain similar results 
including a larger number of vibrational quanta. The excited-state part of the total Hamiltonian, Eq. reads 
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where E is the energy of the optical transition (E = 12210 cm" 1 ), V vv i — J {FC) v {FC) v i , J is the coupling between the 
two BChls and (FC) U = (<?o l e ^) is the Franck-Condon factor characterizing the overlap of the vibrational wavefunction 
in the electronic ground state with that of the v th vibrational level in the electronic excited state. 

In this work, we study the non-rephasing pathways (R\ and R4), and more specifically, their contributions to 
the signal on the diagonal of the 2D spectrum. In particular, we focus on R^ a ' 90 ^ and R^ a ' 91 ^ with a = 2 for 
ujo < 135 cm -1 and a = 3 otherwise (states numbered with increasing energy), because, as will be shown below, 
they both involve the zero- and one-phonon levels local to molecule 1, either on its electronic excited- or ground- 
state manifold. Comparison of the results for these two particular pathways will enable us to discuss the different 
characteristic of the vibronic and the ground-state vibrational coherences, thus taking into account both recent theories 
of the origin of the long-lived oscillations in the 2D spectra of photosynthetic systems [24], [2(| . 

Unlike in Ref. [26], we concentrate on the oscillations of the diagonal peaks in the 2D spectrum. This allows us 
to stay in one-particle approximation [38l |. because, for the ground-state contributions, we are looking at pathway 
jj>(aa,g„) ^ gee [3J3). Since the state \a) is initially excited on the monomer 1 and originates from the ground state of 
the dimer, it cannot contain an excitation of ground-state vibration on the other monomer. Unlike for the crosspeak, 
the one-particle approximation holds on the diagonal cut for all states \a) accessible from the dimer ground state. 



A. Influence of coupling and energy gap on the dipole pre- factor 

As detailed in Section [Vj the initial amplitude of the signal is directly proportional to the dipole pre-factor (see Eq. 
IT4")) . Here, we study the influence of the energy gap AE and site coupling J on the dipole pre-factor, with detailed 
results for pathways R\ involving the zero-phonon level of the excited state on molecule 1 (i.e. eigenstate 1, mainly 
localized on |eigo))- Fig- [3]presents the relative transition dipole moment pre-factor for the dimer response function 
^>(i«.9o) j nvo i vm g electronically excited eigenstate \a) — |2), |3) for ljo = 117 cm' 1 , respective to the pre-factor 
corresponding to the vibrational oscillations on a monomer: 

^ = (4,1 /WSJ • ( 19 ) 

\ / dimer \ / monomer 

Fig. [3^ shows results for the coherence involving the two lowest excited eigenstates |1) and [2). Compared to the 
monomer, the amplitude of the dipole pre-factor in the heterodimer is clearly enhanced (up to more than 8 times) for 
an energy gap comparable to the vibrational energy. It should be noted that, in the range of very small coupling the 
excitonic model breaks down, and the electron-phonon coupling might effectively destroy the excitonic mixing. For 
moderate resonance coupling values and increasing its value, the region for which the amplitude is enhanced through 
borrowing of transition dipole moment spreads over a wider range of energy gaps corresponding to a wider mixing 
region. Note that the studied excited-state coherence has a prevailingly vibrational character only in the shaded area 
of Fig. Hi. 

Fig. (2b presents the amplitude of the dipole pre-factor in the coupled dimer relative to the monomeric one for 
coherences p\ a as a function of the energy gap AE, using Jo = —53.5 cm -1 . We also present the local vibrational 
character \ °f the coherences (c/ Eq. IT5|) . We can define a characteristic energy gap AE X — 94 cm -1 for which 
the coherences pi2 and P13 are equally delocalized, i.e. X12 — Xi3- The highest enhancement (~ 6.5 times higher 
than that in a monomer) is obtained for the coherence involving the two lowest excitons (1,2). This coherence is of 
prevailingly vibrational character for AE > 99 cm -1 . The coherence pis exhibits only a low enhancement (< 2) in 
the region where is it prevailingly vibrational (10 < AE < 88 cm" 1 ) and will therefore make minor contributions to 
long-lived coherences in the 2D spectrum. 



B. Amplitude of the 2D signal and life time of coherences 

In this section, we study the response functions R\ and R4 involving the vibrational coherence in the electronic 
excited-state and in the electronic ground-state manifolds, respectively. The calculations have been performed on the 
FMO inspired dimer at 77 K for various vibrational frequencies, using an overdamped Brownian oscillator mode bath 
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Figure 3: Relative amplitude of the dimer dipole moment pre-factor over that of the monomer (Ar\ at g Q ) for the non-rephasing 
pathway Ri at initial time (a) involving coherence (1,2) as a function of the energy gap AE and coupling J, (b) involving 
coherence (l,a), with a — 2 (red) and a — 3 (magenta), as a function of the energy gap AE for cjq = 117 cm -1 . The plain 
thick lines represent the amplitudes, the dashed lines denote the associated measure of the vibronic character (0.1 Xic*)- In 
both sub-figures, the prevailingly vibrational character of the eigenstates (xia > 50%) is highlighted with colored areas. A 
significant enhancement of the dipole moment pre-factor can be observed for coherence involving eigenstates 1 and 2, which is 
strongly vibrational for AE > 99 cm - . 



(see e.g. (35|) with the Debye frequency A = 130 cm -1 and reorganization energy A = 35 cm -1 . Pigment energies 
were sampled from a Gaussian distribution with 80 cm -1 FWHM centered at the reference energy gap (AEq = 110 
cm -1 ). The energy transfer rates that contribute to the dephasing of the coherences were calculated from the same 
energy gap correlation functions as the line-s hap e functions (see Appendix) using the standard Redfield theory. The 
theory is identical with the one used in Ref. |24j . 

In Fig. [5J we present the energy gap dependence of the diagonal cut (u>i — cl> 3 ) through the Fourier transform of 
the response functions R\ and i?4 before accounting for the energy disorder 

S$ go (u u AE) = A^JAE) G lgo (ui,AE) max (G' 12 (t 2 , AE)) G lgo (u u AE), (20) 
S^^AE) = A^ gi (AE) G 2go (u>i, AE) max (G' gigo (t 2 , AE)) G 2go (wi, AE). (21) 

The energy of eigenstates 1 and 2 depends on both AE and luq (ljq is fixed at 117 cm -1 here). There is a clear 
resonance occurring when the energy gap AE is comparable to the frequency of the nuclear mode, AE ~ ujq. This 
resonance demonstrates the borrowing of dipole moment from the electronic to the vibrational transition, as suggested 
by [H| and illustrated in Section |TTJ The oscillating signal amplitude also depends on the line-shape function G(uS) 
which now also depends on the mixing and therefore on AE. This dependence leads to an even larger enhancement 
of the oscillation amplitude (with respect to the one observed for the vibrational oscillations on a monomer) than the 
transition dipole moment pre-factor alone would predict. 

The complete theory of enhancement has to include averaging over the distribution of energy gaps to account for 
the energy disorder. Fig. [5] presents the maximal amplitude of the signal from pathway R\ involving coherence p\ a 
relative to that of the purely vibrational coherence in a monomer (A^\ defined below), as well as that from R4 (A^), 

m } = max ;5" o(u;i / 4tVl)]dimCT , (22) 

ret 1 ) ( .(max) \i > ' 

maX Pl2,sov' U l' t 2 ' w l ) I monomer 

ro( 4 ) f .(max) 11 
.(4)/ \ _ maxp aa , gi ^i,t 2 ,^ljjdimcr ,„„s 

T ~ m QV [Q( 4 ) (, , -f( max ) , ,„Y| ' 

maxp 22 gi (Wi, t 2 , U)\ JJmonomor 

along the diagonal cut of the 2D spectrum (013 = u)\). The letter a denotes the exciton level which is composed 
mainly of the first vibrationally excited level located on molecule 1 (a = 2 for luq < 135 cm -1 and 3 otherwise), tf* ' 



represents the time of the first maximum after a time t' (t 1 — or 1 ps in the following) and i^ nax ^> j s the time at which 
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Figure 4: Absolute value of the oscillating signal amplitudes S{ 2 go an d S22 gi (Eqs. (I20|l and (J2TJ ) for pathways Ri (a) and R4, 
(b) as a function of the energy gap AE and excitation frequency wi . In both pathways, the amplitude is significantly enhanced 
in the vicinity of the vibrational mode energy, i.e. for AE ~ wo (wo = 117 cm -1 here) due to the resonance occurring. 




-100 -50 50 100 40 60 80 100 120 140 160 180 200 40 60 80 100 120 140 160 180 200 

(Dj-G^p,. (cm 1 ) Nuclear mode frequency O) (cm" 1 ) Nuclear mode frequency co (cm" 1 ) 

Figure 5: (a) Relative amplitude of the signal A r at initial time for pathways _Ri and R4 (Eqs. (|22[) and (|23[) ) involving the 
prevailingly vibrational coherence pi2 as a function of the excitation frequency wi for wo = 117 cm -1 . The enhancement of the 
vibronic coherence (pathway R\) is up to 3 times larger than that of the ground-state vibrational coherence (pathway R4). (b, 
c) Relative amplitude of the diagonal signal through pathways Ri (b) and R4, (c) at wi = w op t as function of the vibrational 
mode frequency wo for a — 2 (red, blue) and a — 3 (magenta, cyan). The shaded areas represent the domains in which the 
coherence is enhanced and of prevailingly vibrational character. The dotted dashed lines show the vibronic coherence life time, 
and the dashed lines represent the signal relative amplitude after 1 ps. The inset in (c) presents the frequency of the oscillating 
signal through the different pathways. The vertical black dashed line at wo — 135 cm - delimits the mode frequency at which 
the ordering of the eigenstates is reorganized and the vibrational character of the coherences changes. It appears that long-lived 
(>1 ps) vibronic coherences are created for wo < 120 cm -1 (b) whereas any mode generates long-lived ground-state vibrational 
coherences (c). 

the signal is maximum for non-decaying coherences. Fig. [5^ presents results for luq — 117 cm^ 1 (a = 2) as a function 
of the excitation frequency lu\ relatively to the frequency w opt , which corresponds to the absorption maximum (for 
luq = 117cm -1 , o; op t = 12165 and 12300 cm -1 for R\ and R4, respectively). For both pathways, the signal amplitude 
in the coupled dimer is clearly enhanced compared to that of the corresponding monomer (maximum 5 times higher 
for i?4, and more than 16 times for R±). Fig. [S£i also shows the life time of coherence P12 involved in obtained 
from a fit of the calculated signal - the life time of the ground-state coherence (pathway R4) is not displayed because 
the purely vibrational coherences are assumed non-decaying as in Refs. [U,[2g. Using ojq = 117 cm -1 , the observed 
excited-state coherence appears to survive for slightly more than 1 ps. 

Fig. [SJd shows the maximal enhancement obtained at initial time (plain lines) and after 1 ps (dashed lines) for 
different vibrational frequencies (on the diagonal cut at uj\ = w op t in the 2D spectrum) for pathway R±. Results for 
pathway i?4 are presented in Fig. [SJ:. In order to account for the dependency of the eigenstate composition on the 
mode frequency and the resulting reordering of eigenstates at loq — 135 cm , results are presented for coherences 
P12 and P13 involved in both pathways R\ and R4. For the relevant domain in which the vibronic coherence 
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Table I: Composition of the eigenstates and eigenenergies averaged over the energy disorder for = 117 cm l . Significant 
contributions are highlighted with bold font. 
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is of prevailingly vibrational character (x > 0.5), is highlighted with shaded areas in Fig. \5jp. For pathway R4, the 
ground-state vibrational coherence remains purely vibrational independently of the origin of the excited state which 
participates in the pathway. For both pathways, the enhancement is larger when coherence pi2 is involved. Comparing 
Figs. [SJd and [SJ;, it can be seen that, at initial time, the enhancement for pathway R\ is significantly (up to 3 times) 
larger than that for R4. Also, for both pathways, there is a broad resonance in the enhancement in the vicinity of the 
energy gap value AEq = 110 cm -1 , which is in accordance with the similarly broad resonance for the enhancement 
of the ground-state contribution for a crosspeak reported in Ref. [2(| . Because the excited-state coherences involve 
mixing of electronic and vibrational states, their life time depends, among others, on the dephasing time of the excited 
states. Consequently, long-lived (>1 ps) coherences will only be created for ujq < 120 cm -1 for the pathway R\. In 
the case of the pathway R4. any nuclear mode will generate long-lived coherences because the purely vibrational 
coherences decay slowly. Therefore for higher mode frequencies, long-lived coherences should originate from pathway 
i?4 only. We studied in detail the R\ pathway involving P12 using ujq = 185 cm -1 (not presented here). We concluded 
that, although at initial time the signal is of comparable amplitude with loq = 117 cm -1 , it is almost absent after 1 ps 
due to the decay. The inset in Fig. [SJ: shows the frequency of the oscillating signals. It is verified that the ground-state 
vibrational coherences oscillate at the frequency of the nuclear mode, whereas the frequency of the vibronic coherences 
does not exactly matches it due to the excitonic splitting effect and the resulting dependence on the site coupling and 
excited-state energies. 

Table U presents the contributions of the local basis excitations to the vibronic states averaged over energy disorder 
for ljo = 117 cm -1 . Eigenstate 1 consists of 85% excitation of the v — transition on site 1 and 15% on site 2. 
Eigenstate 2 corresponds to 64% of the v = 1 transition on site 1. This composition explains the results presented 
above, and confirms that prevailingly vibrational coherences exhibit a prolonged life time. 

Figs. [5^ and [SJd show the evolution of the signal (at uji — u opt ) resulting from pathway R\ and involving the 
prevailingly vibrational coherence for two different nuclear mode frequencies (loq = 117 and 185 cm" 1 ), with different 
energy gaps, namely the reference energy gap (AEq = 110 cm -1 ), the resonant condition AE — u>q, and AE sampled 
from the Gaussian distribution centered on the reference energy gap (denoted by (. . . ) A in Fig. [5]). The fitted 
frequencies and life time of the signal are indicated in the figures, along with the vibrational character \ °f the 
coherence. For the mode frequency luq = 117 cm -1 , we can expect long-lived oscillations with large amplitudes for 
any of the presented energy gaps, because the latter is either equal to or neighboring the vibrational mode frequency. 
Fig. [5^ confirms this expectation. A different behavior is observed for ljq = 185 cm -1 . Here, the amplitude of 
the vibronic coherence is significantly increased (from 0.4 to more than 10 times that in a monomer) when the 
energy gap corresponds to the frequency of the nuclear mode, which illustrates the resonance interaction between 
the two monomers. However, because the averaged energy is sampled around 110 cm -1 , the total signal only slightly 
benefits from the enhancement mechanism and its amplitude is therefore smaller compared to that obtained using 
ujq = 117cm . 

To sum up, the study of a model dimer has enabled us to identify the resonance between the electronic energy gap 
and the vibrational frequency as a crucial condition of the observation of long-lived oscillations in 2D spectra. The 
enhancement mechanism lies in intensity borrowing by excitonic states with strong vibrational character (the vibronic 
excitons) from the strongly allowed electronic states. The vibronic states thus exhibit a significant intensity despite 
the small Huang- Rhys factor of their vibrational component. This mechanism leads to an enhancement of the initial 
amplitude of both vibronic and ground-state vibrational coherences over a bandwidth of about 100 cm - 1 . However, the 
mixing of the electronic and vibrational character of the vibronic coherences implies that these coherences experience 
additional dephasing. The strong dependence of the life-time enhancement on the mode frequency means that vibronic 
coherences will be observable at long population times only when the frequency of the vibrational mode is close to 
the electronic resonance (cf Figs. [SJd and[S]). For the present model, this dephasing leads to similar amplitude of 
the vibronic and ground-state vibrational coherences on a picosecond timescale despite the stronger enhancement of 
the initial amplitude of the former. It does not therefore seem to be possible to exclude one or the other type of 



12 



(a) d) =117 cm' 1 (b) Q) =185 cm" 




012340 1 2 3 4 

Time t 2 (ps) Time t 2 (ps) 



Figure 6: Time evolution of the signal involving coherences of prevailingly vibrational character originating in pathway Ri for 
two mode frequencies: uio — 117 cm -1 (a) and 185 cm" 1 (b). The dynamics of the signal has been computed with different 
energy gaps, namely the energy gap averaged over the energy disorder (denoted by (. . . } A ), the resonant condition (AE — cj a ) 
and the reference gap energy (AE = AEo = 110 cm -1 ). The signal amplitude is presented relatively to that of the monomer 
in the region where it is maximum (wi = w pt)- In part (a) the mode frequency in the vicinity of the energy gap (ujq = 117 
cm -1 ) creates long-lived oscillations of significant amplitude even with energy disorder. In part (b), i. e. for a mode frequency 
out of resonance with the energy gap (too = 185 cm" 1 ), the signal benefits only sightly from the enhancement mechanism, and 
it remains of low amplitude. 

coherence a priori based on amplitude arguments. Both mechanisms will be most likely observed simultaneously in an 
experiment. The most useful difference between the two types of coherences is the difference in oscillation frequency 
(cf inset of Fig. [5}:). The vibrational coherence oscillates with a frequency which is equal to that of the vibrational 
mode, while the frequency of the vibronic-exciton coherence is shifted due to the excitonic interaction. 

VII. DISCUSSION 

The detailed investigation of the interaction between nuclear and electronic DOF in a model dimer provides more 
insight into the mechanism of enhancement of the amplitude and life time of the oscillations seen in 2D experiments. 
These results serve as a verification of the mechanism proposed by Christensson et al. (24|, and form a basis for the 
discussion of this model in relation to the one proposed by Tiwari et al. (26|. We point out that despite technical 
differences, such as the subtraction of the symmetric linear combination of the nuclear mode coordinates in Tiwari et 
al. [26j . the Hamiltonian used in these two works is the same, and so is the basics of the enhancement mechanism. 
The contribution of vibrational states to the optical signal is enhanced by their interaction with electronic DOF. The 
models differ only in the manifold in which the oscillations take place. Christensson et al. [24| suggested the excited 
state, while Tiwari et al. (2(| studied the ground-state contribution. As we have illustrated above, the amplitude of 
both ground-state vibrational and vibronic coherences is significantly enhanced in a system of coupled molecules as 
compared to the isolated monomers. The recent 2D experiments on BChls in solution did not find any significant 
vibrational coherences |39j. These results are in line with what one would expect based on the Huang- Rhys factors of 
the low frequency vibrational modes determined by previous experiments. However, in aggregates, an enhancement 
by factor of five at 1 ps was found for both excited- and ground-state contributions, strong enough to elevate the 
signal above the noise level. 

The ability to distinguish ground- and excited-state contributions to non- linear optical signals is notoriously difficult, 
and to identify unique signatures of electronic, vibronic-exciton and vibrational coherence is less than trivial (40l — 
|43| . One promising proposal has been the use of specific polarization sequences. It has been shown that a specific 
combination of the polarizations in the four wave mixing sequence can single out coherences involving transitions 
with non-zero angle in the molecular frame jiil |45| . When the vibrational modes are treated as members of the bath 
(i.e. standard exciton picture), vibrational transitions on a single pigment are always parallel, and the orientational 
average for pathways representing vibrational coherences will be zero. Therefore such polarization sequences are 
termed "coherence specific" sequences and have been assumed to only report on electronic coherences [45]. However, 
when the vibrational modes are treated explicitly, the mixing between vibrational and electronic DOF leads to a 
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non-zero angle between the vibronic states, even those largely located on the same pigment. The vibronic coherences 
(pathway i?i) discussed in Ref. [24] will thus not be eliminated by the coherence specific polarization sequence. The 
same argument applies to the ground-state vibrational coherences. When modes are treated explicitly, there will be 
a nonzero angle between /J, ago and fi agi , and pathways of type R4 will also contribute to the signal [26] . We can thus 
conclude that the "coherence specific" sequence is not able to unambiguously distinguish between electronic, vibronic 
and ground-state vibrational coherence. 

To be able to determine the origin of the oscillations observed in the experiments, additional arguments, such 
as the oscillation frequencies or dephasing times, are needed. The ground-state vibrational coherences will oscillate 
with a frequency equal to that of the relevant vibrational mode, while for the vibronic coherences, the oscillation 
frequency will depend on the vibrational frequency as well as on the electronic coupling and site energies. The 
vibronic coherences also experience additional dephasing due to the mixing of electronic and vibrational DOF, while 
the ground-state vibrational coherences only experience vibrational dephasing. These quantitative differences could 
in principle be used to distinguish the two models, if the vibrational frequencies, Huang-Rhys factors and dephasing 
times of the vibrational modes were known. 

The low frequency vibrational spectrum of BChl-a in solutions or in different protein environments have been 
studied by a number of techniques. Studies of isolated BChl-a in frozen solution have demonstrated that the strongest 
vibrational mode is found around 180-190 cm -1 (46T - |49| . These experiments have identified a mode around 160 
cnr 1 , but its Huang- Rhys factor was found to be significantly weaker than the one of the modes around 180-190 
cnr 1 0, S3- Although the results of the experiments are rather consistent, it should be pointed out that the low- 
frequency spectral range is sensitive to the local conformation of the macrocycle and the coordination of the central 
Mg atom. This effect was illustrated by Ratsep et al. [46[ by measuring the FLN spectra of BChl-a in different 
solvents. These experiments demonstrate that results obtained on the isolated pigment might not be a good indicator 
of the frequencies and Huang- Rhys factors observed in the protein environment. 

A widely studies case of "isolated" BChl-a in a protein environment is the accessory BChl in the bacterial reaction 
center. Resonance Raman (RR) spectra of the accessory BChl (B-band) have been reported by several groups [50l — 
|52| . The room temperature spectrum identified modes at 84, 117, 180 and 210 cm -1 , where the latter two were 
considerably stronger [5l|. Measurements at 95 K revealed an additional mode at 160 cm -1 which was ~10 times 
weaker than the mode at 180 cm -1 [5l|. Similar results were obtained at even lower temperature, showing that modes 
around 180 and 220 cnr 1 dominates the spectrum |52| . 

There is only a limited number of studies on the vibrational spectrum of the pigments in the FMO complex. A hole 
burning (HB) study found two modes at 120 and 160 cm" 1 with a comparable strength, but were not able to observe 
any modes with higher frequencies due to the narrow absorption band in FMO [53] . Fluorescence line narrowing 
(FLN), on the other hand, revealed a spectrum which was similar to the RR spectrum of the B-band in the reaction 
center [54l. [55|. The two FLN measurements on the different species are similar, showing a mode at 117 cm -1 and 
several modes between 167-202 cm -1 . The strongest modes in both cases were found around 190 cm -1 . 

Taking all of the data on the low frequency vibrational spectrum of BChl-a, we can conclude that none of the 
reported spectra show the presence of a single strong mode in the frequency range corresponding to that observed 
in the 2D experiments (160 cm" 1 ). Rather, the spectra reveal a large number of vibrational modes with comparable 
Huang- Rhys factors in the range from 80-240 cm -1 corresponding to typical energy splittings in FMO [28]. All these 
modes will thus, to a certain extend, experience enhancement due to the interaction between electronic and vibrational 
DOF. For the ground-state wave-packets, the frequency of the oscillation will match that of the vibrational mode. 
Based on the discussion above, we would predict that the mode at 117 cm -1 experiences the strongest enhancement. 
A rather broad amplitude enhancement curve furthermore implies that multiple oscillation frequencies should be 
observed. For the vibronic coherences, the observed frequency depends on the resonance coupling as well as on the 
vibrational frequencies, the number of pigments and probably also on the number of vibrational modes included in 
the Hamiltonian. For the dimer used here, we predict an oscillating frequency of 132 cm -1 and for the full model, 
the value of 140 cm -1 was found (24|. Increasing the resonance coupling would blue-shift the oscillation frequency 
further and increase the amplitude of the vibronic coherence. However, due to the increased electronic character, these 
coherences would also experience a stronger dephasing and a correspondingly shorter life time. Because the amplitude 
enhancement, oscillation frequencies and dephasing dynamics are all closely related to the electronic structure, a 
simulation including all pigments is required to determine whether it is possible to find a set of couplings which 
results in a strong enhancement, correct oscillation frequency and a long life time. Such a simulation could also 
determine whether vibronic or ground-state vibrational coherences dominate in the case of FMO. Furthermore, it 
would be possible to elucidate why a single oscillation dominates in the 2D spectrum, although multiple vibrational 
modes are expected to experience significant amplitude enhancement according to the reduced dimer models. 

Finally, an enhancement of exciton sizes and energy transport rates resulting from vibronic coupling has recently 
been demonstrated in the interacting pigments of two nearly identical light-harvesting proteins [25j. It was shown 
that the vibronic-exciton model provides a realistic explanation for an enhancement in the kinetics of the system with 
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respect to the purely excitonic model. A robust enhancement of the transfer rates and exciton sizes was observed at 
the resonant condition AE = Huj. The mixing of vibrational and electronic DOF in molecular aggregates is therefore 
implied not only in the enhancement of the observed coherence lifetime, but can also lead to enhanced rates of energy 
transport. 



VIII. CONCLUSIONS 



In this work, we have demonstrated on a model dimer how the interaction between electronic and nuclear degrees 
of freedom leads to enhanced amplitudes of vibronic and ground-state vibrational coherences in 2D spectra. The 
vibronic-exciton model used here also provides, without contradiction with molecular dynamics simulations, a plausi- 
ble explanation for correlated fluctuations in different spectral regions, often postulated in purely excitonic model to 
reproduce long-lived coherences. The enhancement of coherence amplitude requires both sufficient electronic coupling 
and a resonance between the vibrational level on one pigment and the electronic transition on the other. Enhancement 
is found to be even more pronounced when the inhomogeneity of the sites is included. The resonance enhancement 
mechanism affects both the excited-state vibronic-exciton coherences as well as the ground-state vibrational coher- 
ences. For the dimer model used in this work, we find that both types of coherences have a similar magnitude at 
long (picosecond) population times. Whether the vibronic coherences dominate in 2D spectra or not depends on the 
excited-state dephasing and population dynamics as well as the influence of inhomogeneous broadening. The results 
obtained here highlight the properties of the enhancement mechanism. To decide which type of coherence dominates 
in FMO complex would require simulations with all seven pigments and a full spectral density including a larger 
spectrum of vibrational modes. 
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Appendix A: Third-Order Response Function 



In this Appendix, we detail the dipole pre-factors and time-dependent response functions for each Liouville pathways 

The transition dipole moment pre-factors read 

A alg v = (((J-g.a ■ ei){fipg v ■ e 3 ){fi goP ■ e 2 ){ti ago ■ (Al) 

A alg, = ((Vg„P ■ ei)(v ag „ ■ e 3 ){npg ■ e 2 ){ii goa ■ e^n (A2) 

A a/3,g„ = {{^g v p ■ e i){^fig ■ e s,){^ agv ■ e 2 ){^g aa - e x )) n (A3) 

A aig„ = ((VgoP ■ e 4)Ws„ • e 3 )(/%, Q • e 2 )(fx ago ■ ei)) n , (A4) 

where (• • • )n denotes the averaging over a random orientation of the pigments. 
The time-dependent response functions read: 

(A5) 
(A6) 
(A7) 
(A8) 

with (. . . ) A representing the averaging over a random distribution of pigment energies. The evolution propagators of 
the optical and excited-state coherences, G(t) and G'(t), respectively, read 

G ag At) = e-^'^e^'e- 7 ™^'), (A9) 



R (l) 


,{ts>h,ti) = 


(a (1) 


G ctg^ 


(*3 )G' afj (t 2 ) G ago (ti ) \ 




,(*3>*2>il) = 


(a (2) 


G Pg» 


(t 3 )G' 0a (t 2 )G goa (h)) A 


p(3) 
n a(S,g„ 


.(£3)^2) tl) = 


\ A aP,g» 


G Pg» 


ih)G' gagv {t 2 )G goa {h)) A 


R {4} 


,(*3)i2,ii) = 


(a {a) 


Gpgo 


(t 3 )Gl go (t 2 )G ago (h)) A 
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and 



Ql — e -i"> a p t g-(r Q +r p) t e -(7 QQ +7^3-27a^) g(t) ( 



(A10) 



where = £ ° is the energy difference between the two eigenstates a and /3, r a is the eigenstate relaxation rate 
computed from the Redfield tensor and g a p(t) = 7 Q( 3 g(t) denotes the line-shape function. The 7 coefficients can be 
expressed via the expansion coefficients c" m , as 



7a/3 J]] ( C n„m„) ( C n„m„) + J]] ( c n„m ) ( C n^,m ) + ( C n„/mo) ( C rL"io) 



(All) 



Using the Debye spectral density, C"{lo) — ttuj 2 J(uj) = 2A ^ 2 " | ^ V2 , where A is the bath reorganization energy and A the 
Debye frequency, the line-shape function reads, in high-temperature approximation, [35] 



,, ,'2Xk B T .A 



-At 



At- 1) 



(A12) 



with ks the Boltzmann constant and T the temperature. 

The contributions of a particular coherence to the signal in a 2D spectrum is obtained after Fourier transformation 
of the evolution propagators during the time intervals t\ and t%: 



S alg v ( UJ 3,h,UJ 1 ) = { A aLg u G <xg„(u3) G' a p(t 2 ) G ago (ux) 



I (2) 



5 'i 3 /3,9„( tJ 3,i2,wi) = (A^^ gv Gp gv (u 3 )G gQgv (t 2 )Gg oa (uJi) 



,(3) 



I (4) 



A 



(A13) 
(A14) 
(A15) 
(A16) 



where G a g v {u) — / + °° dte 1 ^ U} "su) t r ° f fa a g{t) an( ^ W3 are ^- ne excitation and probing frequencies, respectively. 
After analytical integration, we obtain 



G agu M 



A Q (Aw) 



AA 



[T (A a (Au;))-r(A Q (Aw) , A )] 



(A17) 







where Acj = uj agv — uj and Ao = 7§ — *x are constants, and A Q (w) = ° A 4 " — A . r(^) is the Euler gamma 



function with complex argument which satisfies T(z) = L t z 1 e t dt and r(ev, z) = J" t a 1 e '(it is the incomplete 
gamma function. 
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